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Abstract 

The exclusive radiative B — > K*^/, B — > K^, B s — > cfrf and A& — > A7 decays are studied in 
a new physics scenario with two universal extra dimensions compactified on a chiral square. The 
computed branching fractions depend on the size R of the extra dimensions, and a comparison 
with the available measurements allows us to put bounds on such a fundamental parameter. From 

the mode B° — > K*°j we obtain the most stringent bound: — > 710 GeV. 



PACS numbers: 12.38.Lg, 12.60.-i, 13.20.Hc 



Among the various new physics (NP) scenarios proposed to extend the standard model 
(SM), those with extra dimensions (ED) are particularly interesting, as they might be able 
to solve some of the problems affecting the SM without invoking the existence of new in- 
teractions, changing the geometry of the space-time. Indeed, although the attempt of T. 
Kaluza and O. Klein, as well as of G. Nordstrom, to unify the electromagnetism with gravity 
introducing one extra dimension was unsuccessful [lj, the idea of the possible existence of 
additional space-like dimensions resulted to be fruitful, and has been pursued for sectors of 
fundamental interactions with the formulation of different settings involving ED. Such mod- 
els can provide a unified framework for gravity and other interactions, hints on the hierarchy 
problem, connections with the string theory, interesting dark matter phenomenology and, 
remarkably, they can be at the origin of observable phenomena at the colliders and at the 
flavour factories [2| . A common feature is that the compactification of the extra dimensions 
implies the existence of Kaluza-Klein (KK) partners of SM fields in the four-dimensional 
description of the higher dimensional theory, together with KK modes without SM corre- 
spondents. Models in which all the SM fields propagate in the extra dimensions are usually 
denoted as universal extra dimension (UED) scenarios. 

A simple scenario is the Appelquist, Cheng and Dobrescu (ACD) model |3J, a minimal 
extension of SM in 4 + 1 dimensions, with the extra dimension compactified to the orbifold 
S l /Z 2 and the fifth coordinate y running from to 2nR, y = and y = irR being fixed 
points of the orbifold. In this model the SM particles correspond to the zero modes of fields 
propagating in the compactified extra dimension, and the choice of the extra dimension 
topology is dictated by the need of having chiral fermion zero modes. In addition to the 
zero modes, towers of KK excitations are predicted to exist, corresponding to the heavier 
modes of the fields in the extra dimension; such fields are imposed to be even under a parity 
transformation in the fifth coordinate P5 : y — > —y. On the other hand, fields which are odd 
under P5 propagate in the extra dimension without zero modes, and correspond to particles 
without SM partners. The masses of KK particles depend on the radius R of the extra 
dimension orbifold, the only new parameter with respect to SM. For example, the masses of 

KK bosonic modes are given by = + — — n = 1,2,..., with m the mass of the zero 

R 1 

mode, so that for small values of R, i.e. at large compactification scales, the KK particles 
decouple from the low energy sector. Another property of the ACD model is the conservation 
of the KK parity (— j being the KK number. This prevents tree level contributions 
of Kaluza-Klein states to low energy (// <C processes, forbidding the production of 

a single KK particle off the interaction of standard particles. As a consequence, precise 
electroweak measurements can be used to derive a lower bound to the compactification 
scale: 1/R> 250 — 300 GeV jl]. Moreover, this leads to the possibility that the lightest KK 
particles, for instance the n = 1 Kaluza-Klein excitations of the photon and of the neutrinos, 
are among the dark matter components [5j E]. KK modes could be produced at colliders, 
or indirectly detected through investigations of loop-induced processes getting contributions 
from KK modes. This possibility has been explored considering Flavour Changing Neutral 
Current (FCNC) transitions in the ACD scenario |9j [10], with a stringent constraint 
on 1/R obtained from the exclusive radiative B — > K*^ decay [8] and from the inclusive 
B ->• X s j mode p]. 

Models with universal extra dimensions may have a larger number of additional dimen- 
sions, d > 2, with further theoretical motivations than for d = 1. For instance, a model with 
two compactified UEDs allows to cancel the global SU(2) anomaly with three generations 
[12] ; the compactification of two universal extra dimensions on the so-called "chiral square" 
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induces the suppression of the proton decay rate even with maximal violation of the baryon 
number at the TeV scale, as a result of the preservation of a discrete symmetry (which is 
the generator of a subgroup of the six dimensional Lorentz group) [13J. This is the model 
proposed in [12] and considered here: we denote it as 6UED. The two extra dimensions 
are flat and compactified on a square of side L: < x 4 ,x 5 < L, where x 4 and x 5 are the 
fifth and sixth extra spatial coordinates. The compactification is performed identifying two 
pairs of adjacent sides of the square: (y, 0) = (0,y) and (y, L) = (L,y), for all y G [0, L], 
which amounts to folding the square along a diagonal. The fields are decomposed in Fourier 
modes, in terms of effective four dimensional fields labeled by two indices (I, k). Hence, the 
KK modes are identified by two KK numbers which determine their mass in four dimensions: 
zero modes correspond to SM fields. 

The values of the fields in the points identified through the folding are related by a 
symmetry transformation. For instance, for a scalar field, the field values may differ by a 
phase. The choice of the folding boundary conditions (and of the constraints on such phases) 
is mostly important in the case of fermions, since a suitable choice allows to obtain chiral 
zero modes, while higher KK modes have masses determined (as for scalars) by the relation: 

J\1 _|_ fc2 

Mi t k = — , where R = L/n is now the compactification radius, with I and k integer 

numbers. The theory has an additional symmetry, the invariance under reflection with 
respect to the center of the square, which distinguishes between the various KK excitations 
of a given particle. A KK mode identified by the pair (I, k) of indices changes sign under 
reflection if l + k is odd, while it remains invariant if l + k is even. As a consequence, stability 
of the lightest KK modes is guaranteed, so that such are good candidates for dark matter 
constituents. 

The model describes SM particles, their KK excitations, and new particles without a SM 
correspondent from fields whose Fourier decomposition does not contain a zero mode. As an 
example, two scalar fields originate from the mixing of the fourth and the fifth component 
of the vector gauge bosons. Then, at each KK level, a linear combination due to a further 
mixing of these two fields with the Higgs field is promoted to a pseudo-Goldstone mode 
reabsorbed to give mass to the spin-1 KK mode; another combination named the "spinless 
adjoint" [H] remains as a physical real scalar field. The consequence is an interesting phe- 
nomenology which manifests itself in the production cross sections and decay modes of other 
KK states, which differ from the d — 1 case: as an example, the (1,1) states with mass y/2/R, 
can be produced in the s— channel (T2J [15l [16]. The phenomenology of KK excitations has 
been investigated in view of collider searches [T7J [18] as well as for the cosmological and dark 
matter implications [19J [20] . 

The new particles may contribute as virtual states to FCNC processes, and actually from 

the inclusive B — y X s ^ decay the bound — > 650 GeV (at 95% C.L.) has been obtained 

R 

|21j . The analysis can be extended to the exclusive b — y 57 induced modes, prompted by the 
example of a single UED in which such processes turned out to be the effective in constraining 
the compactification scale. Here we focus on the transitions B —y K*j, B —y 1^1(1430)7, 
B s <Pl and A& — > A7. For some of these processes, precise experimental data are available 
and can be exploited to put bounds on 1/R. The modes that have not been observed, yet, 
can be explored at the CERN LHC and at the new planned flavour factories. 

In SM the effective weak Hamiltonian describing the b — y sj and b — y s gluon transition 



3 



can be written as 1221: 



H, 



G f 



eff 



-^v tb v* [ y £c i o i + c 77 o 77 + c SG o 8G 



(1) 



1=1 



Gf is the Fermi constant and Vij are elements of the Cabibbo-Kobayashi-Maskawa (CKM) 
mixing matrix; terms proportional to V u f,V* s are neglected in (1) since the ratio 



v tb v* s 



is 



O(10 2 ). Cj are Wilson coefficients, and O, are local operators written in terms of quark 
and gluon fields: 
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(2) 



a, /3 are color indices, b R l 



16vr 2 
1±75 



A" 



b, and = -[7^, 7"]; e and are the electromagnetic 



2 2' 

and the strong coupling constant, respectively, m;, and m s the beauty and the strange quark 
masses, F M „ in 77 and G° u in 0§ 7 denote the electromagnetic and the gluonic field strength 
tensors, and A a are the Gell-Mann matrices. 

The Wilson coefficients in have been computed at NNLO in SM |23j. The most 
relevant contribution to b — > S7 comes from the operator 7ll a magnetic penguin specific 
of such a transition. The term proportional to m s contributes much less than the one 
proportional to m&, the reason for which the emission of left-handed photons dominates over 
that of right handed ones in SM. Since the coefficient C 7l depends on the regularization 
scheme, it is convenient to consider, at leading order, a combination that is regularization 
scheme independent [24] : 



where r\ 



at s (pw) 



(3) 



and C^^w) = 1 (the superscript (0) stays for leading log approxima- 



tion), with 



14 



23 ' 23 



16 6 

a 3 



12 



04 = — - 



23 23 

a 5 = 0.4086 , a 6 = -0.4230 , a 7 = -0.8994 , a 8 = 0.1456 , 
ht = 2.2996 , h 2 = -1.0880 , h 3 = ~ , h 4 = , (4) 
h s = -0.6494 , h e = -0.0380 , h 7 = -0.0185 , h 8 = -0.0057 . 
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The ACD and the 6UED models belong to the class of minimal flavour violation mod- 
els, therefore the only modification with respect to the SM consists in a different value of 
the Wilson coefficients in the effective weak Hamiltonian (jl|, without new operators. The 
explicit expression of C$f in the case of two extra dimensions can be found in Ref. [2"T] . 
It should only be mentioned that the sums over the KK modes entering in the expression 
of the Wilson coefficients in the extra dimensional framework diverge logarithmically, and 
should be cut in correspondence of some values of Nkk — I + k, viewing this theory as an 
effective one valid up to a some higher scale. Following [2T], the condition Nkk — 10 can 
be chosen. 

To consider the contribution of the effective weak vertex 0-j 1 to the transitions of interest 
here we need the hadronic matrix elements 

<V{p',r ] )\sa tiV q v b\B {s) {p)> = ie^p rf v p a [p'f 2 if {s) ~ 



(5) 



<V(p',r ] )\sa^ l5 b\B {s) (p)> = rf^{M% — M v ) — {rf ■ q)(p + p') 



+ (v* ■ q) 



M 



B, 



O) 



My 



(P + P) 



rpB{s)- 

1 2 



rpB(s)- 

1 3 



if) 



A' 



(q 2 ) • (6) 



V stays for K* and K\ (1430) in the case of B decays, and for 0(1020) in the case of B s decays; 
q = p — p' is the photon momentum and rj the V meson polarization vector. In the case of 
fQ(1430), which is a spin 2 particle, the polarization vector is described by a two indices 



symmetric and traceless tensor, therefore in (5 



6) r] a = rj 



,0 n 

M f 



Relations exist among the 



three form factors Tf^ v , i = 1,2,3, in particular the condition: T x (s) (0) = T 2 (s) (0). 
Due to this relation, the rate of the processes B^ s ){p) — > V(p', rf)^{q,t) (e the photon 

polarization vector) can be expressed in terms of a single hadronic parameter if (s) ^ y (0) for 
each channel: 
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B x & 
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(8) 



B, 



eff 



167T 2 



and the first equation applies both to B — > K* r y and 



In the numerical analysis we set the particle masses and lifetimes, as well as 
the CKM matrix elements to the PDG values; for the quark masses we use rrib — 4.8 GeV 
and m s ~ 0.130 GeV [25]. For the form factors we use results obtained by QCD sum 
rules [25]: in particular, for B — > K* we use the results obtained by three-point QCD sum 
rules [26], based on the short-distance expansion, which provide Tf^ K * (0) = 0.38 ± 0.06, 
and light-cone QCD sum rules (LCSR) [27], based on the light-cone expansion, which give 
Tf^ K *(0) = 0.333 ± 0.028. These values are larger than the lattice QCD result obtained 
in quenched approximation [28]. LCSR calculations of B — > X|(1430) and B s — > (j) form 
factors give: if ^ A ' 2 *(0) = 0.17 ± 0.03 ± 0.04 [29] and if a ^(0) = 0.349 ± 0.033 ± 0.04 [27]. 
In the case of — > A7 we define the matrix elements 

< A(p',s')|w^6|A b (p,s) >=u K [f^{q 2 ) lil + if^{q 2 )a^q u + fl{q 2 ) qil ]u K (9) 
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Table I: PDG averages for the branching fractions of radiative B, B s and A& decay modes |34j . 



mode 


BR 


B + K *+ 1 


(42.1 ± 1.8) x 10" 6 


B° -> K*° 7 


(43.3 ± 1.5) x 10~ 6 


B + -> K|(1430)+7 


(14 ±4) x 1(T 6 


B° -> Iff (1430)° 7 


(12.4 ±2.4) x 10 -6 


B s -> 07 


(57±?g) x 1(T 6 


A fe -> A 7 


< 1.3 x 1(T 3 (90% C.L.) 



< Hp', s')\sia^q u -i b b\X b (p, s) >= U A {g({q 2 )^^ b + ig 2 (q 2 )a^ u q"-f 5 + gj(q 2 )q^ 5 ] u Ab (10) 
with «a and u\ b the A and A& spinors; s denotes the baryon spin. The determinations of 



the form factors in Q9j)-( 10 ) are quite uncertain. However, it is possible to invoke heavy 
quark symmetries for the hadronic matrix elements between an initial spin=^ heavy baryon 
comprising a single heavy quark Q and a final spin=| light baryon; due to the heavy quark 
symmetries the number of independent form factors is two, since for tuq — > oo and a generic 



Dirac matrix F one can write [30] 

<A(p',s')\sTb\A b (p,s) >=u A (p , ,s'){F 1 (p' -v)+ jt F 2 {p' ■ v)}Tu Ab {v, s) (11) 
where v = jjj^- is the A b four-velocity. The form factors F\ )2 depend on the invariant 



M\ +M\-q 2 

p ■ v = — — , but we refer to their q dependence for convenience. Using Eq.(jll|) 

-•> 'A, 

we can relate the form factors in Eqs. ([£])-( 10) to the functions Fx^m Eq.(ll), 



fi = 

^ = -( 1 -^T L ) / -2 ^ 



T 

9 3 = 

at momentum transfer close to the maximum value q 2 ~ q%^ ax = {M\ b — Ma) 2 . We assume 
their validity to the whole phase space, introducing a model dependence in the predictions. 
The Ab — > A 7 decay width reads in terms of the form factors F± and F 2 : 




r(A b ^A 7 ) = ^(FU0) + F 2 (0)^ K + m 2 )M| b (l-^) . (13) 

A determination of F\ and F 2 has been obtained by three-point QCD sum rules in the 
rriQ — y oo limit [31J. In the following we use F\, F 2 worked out in [10] updating some of the 
parameters used in [3]]: Fi(0) = 0.322 ± 0.015, and F 2 (0) = -0.054 ± 0.020 [10]. Other 
determinations of the form factors in the transition A b — > A 7 have been performed in [32] 
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Figure 1: Predicted branching fractions of B + — > K* + j (upper panels) and B° — > K* ^ (lower 
panels) as a function of the compactification parameter 1/R (in units of GeV) , using the form factor 
in |26j (left) and in [27] (right). The horizontal bands correspond to SM theoretical expectations 
with 1 a uncertainties (yellow [light]) and experimental measurements with 2 a uncertainties (blue 
[dark]) respectively. 



at large recoil and in [33] in the covariant constituent quark model. We set the and A 
masses and the lifetime to their PDG values [31]. 

The computed branching fractions are functions of — , through the dependence of the 

coefficient C-j in the 6UED model. The results are depicted in figsJTJ |2J where we have 
included the uncertainty on the form factor value at q 2 = 0, and a second uncertainty, 
intrinsic of the model, due to the choice of the matching scale in the calculation of C-j and 
of the value of Nkk- These last uncertainty is discussed in [21], together with another 
one which comes from fixing two boundary couplings hip which are 0(1) and enter in the 
expression of the masses of the Higgs fields in this model. Altogether these uncertainties do 
not exceed |21j . and we include this range in our error on the branching fractions. The 
experimental data for the various branching ratios are collected in Table [I] and represent 
PDG averages [31]. For B — > K*^ modes, the results represent averages of BaBar [35], Belle 
[36] and CLEO [37] measurements. For B — > the result is determined on the basis 
of the analysis in [38], while for B s — > (f) it stems from ref. [39] . The upper bound on the 
23 (At, — > A7) has been obtained in [H]. Experimental data are represented as horizontal blue 
[dark] bands (at 95% c.l.) in fig{TJ and in fig|2]for B° — > K^ ^ (due the large experimental 
uncertainty, we only show the neutral mode where the error is smaller) and B s — > (fry. 

In the two upper plots in fig{T]S(i? + — > K* +/ y) is computed using either the form factor 
Ti in [26] (left panel) or that derived in [27] (right panel) . The same applies to the two lower 
plots, where the neutral channel B° — > K* 0/ -y is considered. There is a model dependence, 
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1 



1 



with the resulting bounds : — > 397 GeV (charged channel, form factors in [26]), — > 564 

R R 

GeV (charged channel, form factors in [27]), — > 433 GeV (neutral channel, form factors 

1 

in [26]), — > 710 GeV (neutral channel, form factors in [27]). As for B° — > K% 0r y, we obtain 
R 

' " " " 'refer to B s and A& decays. In a previous analysis 



R 



> 324 GeV. The other two plots in fig 

i(0 



0] also the decays B — > Kif'^ has been considered in order to constrain UEDs models 

(with 1 and 2 extra dimensions respectively), finding for 6UED a lower bound of — > 400 

R 

GeV. 




600 
l/R(GeV) 



1000 



Figure 2: Predicted branching fractions of B — > K2J, B s — > (jyy and Aj, — > A7 as a function of the 
compactification parameter 1/R (in units of GeV). The SM results (with 1 a uncertainties) and 
experimental measurements (with 2a uncertainties) are represented as horizontal yellow [light] and 
blue [dark] bands, respectively. 



Bounds from direct searches of KK modes have been discussed in [T7], where the hadron 
collider phenomenology of (1,0) KK modes in the 6UED model was studied. The limit 

— > 270 GeV was found, set from direct searches at the Tevatron. Our bound is more 
R 

restrictive than the one obtained from the inclusive radiative B transition [21], as in the 
case of a single UED. Since the discovery signals of UEDs at the CERN LHC or at a 
future e + e~ linear collider, searching for leptons plus missing energy, are expected with the 
experimental reach of about 1/R ~ 1 TeV [HI H21 S3] , the complementarity of the studies 
of indirect effects of possible UEDs in flavour observables is noticeable. 
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